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ABSTRACT 

U. S. Navy ships not equipped with NTDS currently per- 
form surface radar target tracking manually, a tedious and 
inaccurate process. 

This thesis establishes that an 8-bit microprocessor, 
interfaced to a common surface search radar, could compute 
the course and speed of multiple radar targets. The 
investigation includes the estimation of radar measurement 
errors, a determination of practical digital interface 
limitations, and the development of a Kalman filter 
algorithm, a floating point arithmetic library for the 
INTEL 8080 microprocessor, aind the microprocessor tracking 
system software. The error statistics and execution time 
of the microprocessor software are presented for several 
computer simulated target tracks. 
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I . INTRODUCTION 



"In December 1973 INTEL Corporation shipped the first 
8-bit N-channel microprocessor, the 8080." /~1_7 In the 
months that followed, several semiconductor manufacturers 
introduced their own version of INTEL'S 8080 microprocessor 
The resulting competition has driven microprocessor prices 
down and performance up. Advanced Micro Devices annoimced 
in April 1976, their AM9080A-4 microprocessor, with an 
instruction execution time one-half that of INTEL'S 
original 8080, was available for $21.00 /~2_7. In the 
same annoimcement , they predicted the AM9080A would cost 
$5*00 by 1980. Paralleling the advances in microprocessor 
technology, the manufacturers of semiconductor memories are 
making steady progress in expanding the memory capability 
of single semiconductor chips, reducing memory access time 
and subsequently reducing the cost of high-speed semiconduc 
tor memories used with microprocessors. 

These advances in semiconductor technology are taking 
place at a time when the U. S. Navy is debating the cost 
of rebuilding the fleet. A significant portion of the cost 
of building a modern warship is the electronics in its 
weapons systems. Microprocessors offer the potential to 
(1) reduce the cost of digital systems, (2) perform complex 
functions at remote stations, relieving the congestion at 
larger central computing facilities, and (3) perform func- 
tions currently handled by watch personnel, thus reducing 
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the manning requirements of watch sections. 



A. OBJECTIVES OF THESIS 

The primary objective of this thesis is to demonstrate 
the capability of microprocessors to perform surface search 
radar tracking, thus relieving the tedium of this function 
from the typical underway watch section. 

The secondary objective of this thesis Is to investigate 
the performance of a Kalman filter in the two-dimensional 
polar coordinate system, with measurement noise, typical 
of a surface search radar. 

B. SCOPE OF THE WORK 

The scope of the work was limited to the development of 
a Kalman filter tracking algorithm for the INTEL 8080 
microprocessor and simulation of circuitry necessary to 
interface with the AN/SPS-10 surface search radar. 

A Kalman filter was chosen because it is computationally 
the most complex of the tracking filters typically used 
for radar tracking. Thus, it would be the most demanding 
on the microprocessor. Additionally, the Kalman filter is 
seldom implemented in its complete form due to its computa- 
tional complexity thus adding additional relevance to this 
investigation. 

The parameters for the AN/SPS-10 surface search radar 
were used for computations requiring specific radar para- 
meters. This radar was selected because it was available 
at the Naval Postgraduate School and is installed aboard 
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the majority of U. S. Naval ships. 

INTEL Corporation's 8080 microprocessor was used for 
the implementation of the Kalman filter because it was 
available at the Naval Postgraduate School along with sup- 
porting software. 

The interface between the microprocessor and the radar 
was not constructed due to: (1) Cost, time and manpower 

constraints, and (2) The performance of the system could 
not be effectively measured against non-cooperative targets 
of opportunity. 

The performance of the radar and the interface equipment 
were effectively simulated on the Naval Postgraduate 
School's IBM -360 computer, thus permitting comparison of 
the microprocessor's performance against that of the IBM- 
360 system for known target tracks. 

The output of the tracking system would be displayed 
for human utilization rather than being eleotrically 
interfaced to fire control equipment. The minimum standards 
of performance were established as that obtainable from 
a well trained, fully manned underway watch section typical 
of that employed on a modem destroyer not equipped with 
a Naval Tactical Data System (NTDS). 
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II . BACKGROUND 



A. PERFORMANCE OF A TYPICAL CIC WATCH TEAM 

Onboard naval ships not equipped with Naval Tactical 
Data Systems (NTDS), the tracking of radar contacts, both 
surface and airborne, is performed manually by the Combat 
Information Center (CIC) watch team. During normal peace- 
time steaming, this watch team may consist of from two to 
ten personnel, depending on the size of the ship and the 
complexity of its CIC equipment. Due to the risk of col- 
lision with other ships, their effort is primarily devoted 
to detecting, tracking and computing the course, speed, 
and closest point of approach (CPA) of other ships within 
a 20 to 30 nautical mile radius. This can be a tedious, 
error-prone task which generally requires a duplicate effort 
by the bridge watch team in order to avoid the awesome 
consequences of a collision at sea. The installation of 
an economical device to track and compute the tactical 
data for surface radar contacts could eliminate one or more 
watch-standers in CIC while improving the reliability of 
the information provided to the bridge watch team. 

The manual method of computing the course, speed and 
CPA of a radar target involves a graphical solution of the 
true and relative velocity vectors of the radar target's 
motion and own ship's course and speed. The target's 
relative position measurements are obtained from a radar 
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repeater. The range measurement accuracy is at "best + 20 
yards with a well calibrated repeater and a well trained 
operator. The bearing measurement accuracy is at best 
+ 1 degree with a well trained operator. The own ship's 
course and speed used in determining the own ship's velocity 
vector are normally based on the ordered course and speed. 
The course actually steered by the helmsman may vary from 
the ordered course by + 2 degrees. Depending on the condi- 
tion of the ship's hull and its engineering plant, the 
speed will be accurate within + 1 knot. The combined 
results of these errors lead to a target course accurate 
to within + 10 degrees and a target speed accurate to 
within + 3 knots, not allowing for human error. The major 
disadvantages of the manual solution, in addition to its 
poor accuracy, is the time required to obtain the solution 
and the inability to obtain solutions if either the target 
or own ship maneuvers between radar position measurements, 
typically 3 or more minutes apart. The estimates of error 
are based on the author's personal observation of unalerted 
CIC watch sections. 

B. THE TRACKING PROBLEM 

The "Tracking Problem" as referred to in this investiga- 
tion is that of using radar measurements to accurately 
determine the course and speed of a waterborne radar target. 
The targets of interest include merchantmen, warships, and 
other ships capable of open ocean operations. As such. 
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their motion on the surface of the ocean can he approximated 
by a linear dynamic system. In a sampled data system, such 
as a search radar where measurements are taken once each 
antenna rotation, the target dynamics can be represented 
in matrix form as Z~3_7* 
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Z_ = Vehicle State at Scan n 
n 

- E/W Postion 

» a 

- E/W Velocity 

- N/S Position 
Y^ - N/S Velocity 



T = Time Between Scan n and n+1 



a = Acceleration Acting on the Target from Scan n 
to n + 1 (E/W and N/S components) 



If the radar measurements were perfectly accurate and 
the target had zero acceleration, then the tracking problem 
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would be a simple velocity vector problem (Figure 1). 
Unfortunately, the radar bearing and range to the target on 
any particular scan is contaminated by errors which are a 
function of the radar system and the received signal to 
noise ratio. The measurement errors, to be discussed in 
the next section, are characteristically Gaussian. An 
optimal estimator is well suited to tracking problems such 
as this . 

"An optimal estimator is a computational algorithm that 
processes measurements to deduce a minimum error estimate 
of the state of a system by utilizing: knowledge of 

system and measurement dynamics, assumed statistics of 
system noises and measurement errors, and initial condition 
information." 

A Kalman filter is a common optimal filtering technique 
for estimating the state of a linear system is 

chosen for implementation because it is computationally the 
most demanding as compared to other techniques used in 
similar tracking systems and, as such, places the most 
demand upon the microprocessor. The details of the Kalman 
filter will be covered in Chapter IV following an analysis 
of the measurements and associated hardware which provide 
its inputs and determine its parameters . 

C . BtEASUREMENTS 

In order to solve the tracking problem, it is necessary 
to make measurements of (1) target range, (2) target bearing, 
(3) own ship's course, (4) own ship's speed, and (5) time 
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FIGURE la GEOGRAPHIC MOTION PLOT 



Y 




16 



Y 




FIGURE Ic VELOCITY VECTOR PLOT 
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of the measurement. It is assumed that all noise is Gaussian 
and received signal to noise ratios are large (» 1). 

The target range is determined hy measuring the time 
required for a radar pulse to travel to and from the 
target. 




R - Range to Target 
c - Speed of Light 

t - Elapsed Time from Pulse Transmission to 
Pulse Detection 

Using threshold detection to obtain a leading edge measure- 
ment, the theoretical standard deviation of error for the 
measurement of t is /~5_7* 

i Tr ■ ( 

= T ^ 

(2 s/n)^ ^ 

At = Error in Elapsed Time Measurement t 

= Rise Time of Received Pulse 
r 

s/n = Signal to Noise Ratio of Video Pulse 

If the rise time of the pulse ( Tj,) is limited by the receiver 
IF bandwidth B, then ^ Z~5_7* The theoretical 

standard deviation of range error using leading edge detec- 
tion is therefore: 

c ^ 

°R " 2B(2 s/n) 2 

cr^ = Standard Deviation of Range Error 
The mean range error is negligible when compared to the 
range measurement. Additionally, the mean error will tend 



( At)^ 



18 



to be cancelled by subtraction in the relative velocity 
calculations . Thus the computed target course and speed 
will not be affected by the mean range error. 

The target bearing is determined by measuring the bearing 
of the radar beam at the time the target video is detected. 
The mean error of the bearing measurement will be a fimction 
of the antenna beam pattern, and the target range and radar 
cross section. The resulting bearing to target will not 
be through the center of the target but the computed target 
course and speed will not be affected. 

Using beam splitting, a technique for estimating the 
center bearing of the target, the theoretical minimum stan- 
dard deviation of bearing error for a Gaussian beam distri- 
bution is Z~5_7* 



where 



<^0 



1.06 0 ^ 

Ng(s/n)^ 



I®"'' U.C) 

V 



<f^ = standard Deviation of Bearing Error 
0 = Two Way Beamwidth (= 0.4250 0 -q = Half 

g D D 

Power Beamwidth) (Same units as Oq) 

(s/n) = Signal to Noise Ratio at Center of Beam 

N = Number of pulses transmitted as the antenna 
° rotates through 2 0 

& 



The standard deviation of both range error and bearing 
error are functions of the video signal to noise ratio. 

The video signal to noise ratio for a radar with single 
pulse detection is given by 
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s 



i.U 



n 



(47T)3 k I 3 r'^ 

where 

Ppj^ = Peak Transmitted Power (watts) 

G = Antenna Gain 
X = Wavelength (meters) 

2 

O' = Radar Cross Section of Target (m ) 

ot = Attenuation Constant of Propagation 
Medium (db/m) 

R = Range to Target (meters) 

k = Boltzmann's Constant (1.23 x 10~^^ joule/deg) 

T^ = Standard Temperature (290° K) 

= Receiver Noise Bandwidth (HZ) 

L = System Losses 
s 

When these relationships are applied to the AN/SPS-10 
(series) Surface Search Radar, with the AS-II 6 I/SPS 
Antenna, the standard deviations for bearing and range can 
be related to target range and cross sectional area as 
shown in Figures 2 and 3* The minimum performance charac- 
ters of this radar system are listed below. 



Peak Power = I 90 kw Z”^_7 
Antenna Gain = 30 db 8_7 
Pulse Repetition Rate = 625 HZ Z~6_7 
Frequency = 5825 MHZ Z~^_7 
Propagation Attenuation = .6 x 10“-^ db/m Z~5_7 
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FIGURE 2 STANDARD DEVIATION OF RANGE ERROR VS. RANGE 
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FIGURE 3 STANDARD DEVIATION OF BEARING ERROR VS. RANGE 




I 






Noise Figure = 22.9 db 
Bandwidth = 5 MHZ 

Noise Bandwidth (I .05 x B) Ci 5.25 MHZ 
System Losses 5 db 
Antenna Rotation Rate = 17 RPM 
Antenna Horizontal Beamwidth = 1.9° 



^ 7_7 

Z-5J7 

/■5J7 

ZT8J7 

Zf8j7 



As proposed earlier, a threshold detector could be used 
to indicate the occurrence of the leading edge of the 
received radar video pulse. The threshold to noise ratio 
of such a detector will establish the ranges at which 
various targets can be tracked and the minimxim signal to 
noise ratio which the tracking system will encounter. 
Leading edge detection is preferred because it eliminates 



errors resulting from variations in the pulse width. 

When considering the threshold to noise ratio of such a 
detection scheme, it is necessary to establish the desired 
mean time between false alarms, and the desired minimum 
probability of detection of the target. 

The mean time between false alarms is a measure of 
how often the tracking system will erroneously recognize 
a noise pulse as a target. The tracking system can not 
function effectively if this occurs frequently. A mean 
time between false alarms of 60 sec has been selected.. 

The mean time between false alarms is a function of the 
signal to noise ratio and the IF bandwidth of the receiver 
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(-^Tp-^) = Threshold to Noise Power Ratio 



IF “ 
FA ^ 



IF Bandwidth 

Mean Time Between False Alams 



The exponential tern in this expression is the inverse of 
the probability of false alarm Z~5_7« 



p = Exp(- It ) 

zYo 



. ,1 



C .'b 



,i9 



Thus 



FA T„^ B 



FA ®IF -2 






The mean time between false alams as seen by the tracking 
system can be increased if the tracking system is disabled 
whenever targets are not anticipated. For the AN/SPS-10, 
the maiximum unambiguous range is 129.6 nm of which only 
30 nm is of interest to the tracking system. Therefore, ^ 
the measurement circuitry could be disabled 77 percent of 
the time. This would multiply the mean time between 
false alarms by 4.32. An additional advantage can be 
gained by using bearing and range gates to further limit 
the opportunity for false alarms. If such bearing and 
range gates isolated 1000 yard squares for each target 
and the tracking system were tracking eight targets, the 
mean time between false alarms would be increased by an 
additional factor of l4l4. The combination of these factors 
leads to a probability of false alarm of 2.036 x 10 ^ and 



i£. 
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a threshold to noise power ratio of 10.8 or 10. 3 db. 

The probability of detection of the target is a measure 
of how often the tracking system will recognize the target 
video when a target is present. It is essential to the 
effectiveness of the tracking system that it detect the 
target video reliably. The probability of detection is a 
complex function of the signal to noise ratio, the charac- 
teristics of the detector, and the probability of false 
alarm requiring numerical techniques Z~5_7* Figure 2,7 
of Reference 5 shows a plot of the probability of detection 
of a sine wave in noise as a function of the signal to 
noise power ratio for various probabilities of false alam. 
This plot shows that a signal to noise ratio of 12.5 db to 
14.0 db is required for a- probability of detection of O .90 
to 0.99» respectively, and a probability of false alarm 
of 2.0 X 10“-^ , If all targets were perfect radar reflec- 
tors, it would be possible to predict the performance of 
the threshold detector and the range at which various 
targets could be tracked. However, the reflected radar 
signals from the targets of interest will fluctuate due to 
cancellation and reinforcement of the waves reflected 
from various portions of the target. These fluctuations 
can be modeled as either Swerling’s Case I or Swerling's 
Case II Z~5_7» In either case, an additional signal to 
noise ratio of 8 db to 17 db will be required to insure 
a probability of detection of 0.9 to 0.99i respectively 

zf5_7. 
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Thus the proposed threshold detection scheme should have 



a threshold to noise power ratio of 10.8 to achieve an 
effective mean time between false alarms of 6o sec. This 
would lead to a worse case minimum average signal to noise 
ratio of 20.5 db to 31*0 db for a probability of detection 
of 0.9 to 0 . 99 1 respectively. 

For a radar measurement system with the parameters 
established above, targets with the following average 
radar cross sections could be tracked with the AN/SPS -10 
at the ranges indicated. 



Average Radar Cross 
Section (m^) 

1 

10 



10 ' 



10 - 



10 



10 - 



Range 

(yards/meters ) 

3,600/3,290 

6,400/5.850 

11,500/10,520 

20,000/18,290 

35,000/32,000 

60,000/54,860 



Based on a threshold to noise ratio of 10 . 8 , the stan- 
dard deviation of the radar bearing and range error (Figures 
2 and 3 ) will be O.0763 degrees and 7*057 yards, respec- 



tively. 

The own ship's course and speed measurements come from 
the ship's gyro and pitometer log, respectively. Accordingly, 
their error statistics are a function of the equipment 
installed on the particular ship. This equipment is fre- 
quently associated with the ship's fire control equipment 
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and, as such, is intended for tracking operations. 

The gyro has a dual influence on the tracking prohlem. 

It not only determines own ship course hut also determines 
the true hearing of the radar aintenna. The mean error in 
radar hearing, and likewise own ship's course, will result 
in an error in the solution. This can he viewed as a 
shift in the reference plane which does not prevent the 
correct prediction of the future relative position of 
targets. As a refinement to the tracking system, it would 
he better to measure the average course and speed made 
good since the last measurement rather than the instan- 
taneous values at the time of the measurement. Such an 
average measurement could he achieved hy continuously 
integrating the variations in the resolved N/S and E/W 
velocities using a time constant equal to one-half the 
period of the antenna rotation rate . This integration 
process would have the desirable side-effect of minimizing 
the variance of error in own ship's course and speed. 

For subsequent calculations, it is assumed that the standard 
deviation of error of own ship's gyro is 0.1 degrees and 
own ship's speed indicator is 0.1 knots. 

The fifth and final measurement, the time of the measure- 
ment, could be made with a very high accuracy, depending on 
the quality of the clock used. The objective is to deter- 
mine the elapsed time between position measurements. The 
antenna rotation rate determines the mean elapsed time 
between measurements on a unique target, and the pulse 
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repetition rate of the radar determines the useful precision 
to which the elapsed time needs to be measured. The mean 
rotation rate for the AS-II61/SPS Antenna is 17 RPM. The 
mean pulse repetition rate for the AN/SPS-10 (series) 

Radar is 637 *5 HZ. Thus, the mean elapsed time between 
measurements is 3*53 sec and the useful precision of the 
time measurement is + 0.78^ ni sec. 

D. MICROPROCESSORS 

A microprocessor is a large scale integrated (LSI) 
semiconductor device or family of devices which contain 
an arithmetic logic unit (ALU) and several working registers 
capable of performing binary or BCD arithmetic and/or 
Boolean logic operations. Microprocessors generally fall 
into one of five categories; (1) single-chip controllers, 

(2) 4- bit family, (3) 8-bit family, (4) l6-bit family, 
and (5) bit slices /~9_7* 

The single chip controllers are LSI devices which 
incorporate all the necessary hardware, including program 
and data memory circuitry, on a single chip. Their data 
processing applications include devices such as calculators, 
gas pumps, cash registers, and scales. They include devices 
with 4-bit and 8-bit ALUs. Typically the 4-bit versions 
perform serial BCD operations. They are limited in their 
speed of program execution by the power of their instruc- 
tion sets and the size of their ALUs. The complexity of 
the operations they can perform is limited by the size of 
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their memories. Some manufacturers are developing single- 
chip versions of their 8-bit family devices which will 
expand the capabilities of the single-chip controllers 
to the point of overlapping the 8-bit families 

The if— bit microprocessor families are similar in 
function to the single-chip controllers except the memory is 
not included on the chip. Therefore, the 4-bit families 
are capable of handling larger programs. The term family 
refers to the set of compatible chips which go together 
to make up the microprocessor. They require greater design 
effort than the single chip controllers because several 
chips, including clocks and l/O interface chips, along 
with the separate memory, are required for even the simplest 
task. Like the 4-bit controllers, the 4-bit microprocessor 
families are limited in the speed of their program execution 
by the power of their instruction sets and the size of 
their ALUs. 

The 8-bit families are presently dominating the micro- 
processor applications /~ 9 _ 7 * Their instruction sets are 
approaching the complexity of those used in minicomputers. 
They typically provide for 8-bit data manipulation along 

with l6-bit memory addressing. " INTEL Corporation's 

8080 family, with its enhanced 8 O 8 OA CPU, Motorola Semi- 
conductor's 6800 family, with its enhanced 6800D CPU, 
and Rockwell's PPS-8 family currently rank one, two, and 
three in popularity among users." Z~9_7 

The strength of the 8 O 8 O lies in the power of its 
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instruction set and the extent of its software support. 

The 8080 has ?2 basic instructions, including conditional 
branching, decimal as well as binary arithmetic, logical, 
register to register, stack control and memory reference 
instructions These instructions can be executed 

in 4 to 18 cycles of its 2 MHZ clock allowing some instruc- 
tions to be executed in 2 psec 1 _ 7 « These instructions 
facilitate l6-bit addition and register to register trans- 
fers. Among its software support are a System Monitor, an 
Assembler, a Text Editor, and a PL/M compiler, cross 
assembler, and simulator. The popularity of PL/M, an 
INTEL derivative of PL/l, and the 8080 have caused other 
manufacturers to model their microprocessors after the 
8080, such as Advanced Micro-Devices 9080A, or to develop 
their own version of PL/M such as PL/M 6800 by Intermetrics 
for the Motorola's 6800. 

The 8-bit devices are limited in performance by their 
8-bit data path to and from memory. This necessitates 
using several memory cycles to read or write address and 
other l6-bit variables. 

The l6-bit microprocessor family, including National's 
PACE, Texas Instruments Inc. 's 9900, and General Instrument 
Corp. 's CPI 600 , attempt to reduce the number of instructions 
and machine cycles required to manipulate memory addresses 
(typically l6-bit variables), and data originating from 
A-D converters or going to D-A converters (typically 10-16 
bits). Their performance approaches that of minicomputers. 
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Western Corp.'s MCP l600 was originally designed for DEC'S 
LSI-11 l6-bit minicomputer and was subsequently distributed 
as an independent device /~9_7* 

The bit-slice architectures, such as Advanced Micro 

Devices Inc.'s AM 2900, INTEL Corp.'s 3002, and Texas 

Instrument Inc.'s SBP0400, are taking on the tough process- 

control and high speed controller jobs previously handled 

by minicomputers. A typical bit-slice device is a 2 or 4 

bit vertical slice of a CPU which is horizontally expandable 

in 2 or 4 bit increments to form a 16 bit, 32 bit or other 

multiple of 2 or 4 bit CPU. They utilize Schottky TTL, 

2 

I L or ECL circuitry and thus are much faster than the 
previously mentioned microprocessors which are primarily 
MOS devices. The bit-slice architecture offers the ulti- 
mate in design flexibility while demanding the most effort 
by the designer. The design effort may include an extensive 
amount of microprogramming prior to development of the 
working program. Additionally, bit-slice devices are the 
most expensive form of the microprocessor. 

In summary, the 8-bit and 16-bit microprocessors are 
the best alternatives to the tracking problem. The single- 
chip controllers and 4-bit families have insufficient memory 
and/or are too slow to handle the Kalman filter algorithm. 
The bit slice architectures have excess capability while 
requiring unnecessary design and hardware cost. As will 
be seen later, the tracking problem can be solved accurately 
with l6-bit resolution in the arithmetic operations. Thus 
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an 8-bit microprocessor programmed for double precision 
arithmetic or a l6-bit microprocessor program for single 
precision arithmetic is suitable. 
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Ill . SYSTEM INTERFACES 



The required interface equipment can he grouped into 
three categories; (1) the radar interface, (2) the own 
ship course/speed interface, sind (3) the operator/user 
interfaces. Although the interface hardware was not con- 
structed, it is essential to the validity of the proposed 
tracking system to verify the feasibility of making the 
required radar and own ship measurements within the pro- 
posed tolerances. 

A. RADAR INTERFACE 

The radar interface serves four functions; (1) to con- 
vert the radar video to rectangular pulses compatible 
with the logic circuitry in the range measurement subsystem, 
(2) to identify the radar video of the target whose position 
is to be measured, (3") to measure and store the bearing, 
range and time of measurement of the designated radar 
target, and to transmit a signal to the own ship 

course/speed and the operator/user interfaces to indicate 
when the radar measurements have occurred. 

The conversion of the radar video to a rectangular pulse 
compatible with the bearing and range measurement circuitry 
could best be accomplished at the threshold detector. The 
threshold detector could consist of one of a variety of LSI 
voltage comparators currently available . The input video 
would be in the 0-2 volt range depending on the setting of 



33 



the radar IF aunplifier gain adjustment. The threshold 
voltage would he a technician adjustment made hy first 
setting the IF amplifier gain to obtain a good video 
display on the radar's master repeater and then adjusting 
the threshold voltage on the comparator to just detect a 
known signal. The critical parameters of the threshold 
detector are a bandwidth greater than that of the IF 
amplifier (5 MHZ) and a high noise immimity. 

Target identification could be accomplished either 
with software or hardware . A software solution would result 
in the measurement and processing of a large number of 
targets which the operator had no interest in and would 
place an unnecessary burden on the microprocessor. Thus 
a hardware technique is preferred. A simple solution to 
this problem is the use of bearing and range gates to 
enable the measurement circuitry at the time a particular 
target video is anticipated. 

The width and time of occurrence of the bearing and 
range gates can be determined in several ways, all of which 
are a function of the accuracy with which the future loca- 
tion of the particular target can be estimated. This 
estimate is a function of (1) the maneuverability of the 
target, (2) the elapsed time between measurements and (3) 
the accuracy of the computed course and speed of the target. 

The majority of targets of interest will be capable of 

turning radii of 100-1000 yards, velocities of 0-40 knots, 

2 

and accelerations of 0-.5 yards/sec . The antenna rotation 
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rate determines the mean time between measurements , which 
in the case of the AS-II 6 I/SPS antenna equates to 3*53 
seconds Z~8_7. The accuracy of the computed course and 
speed of the target will depend on the tracking system 
parameters and the maneuverability of the target. 

After the three major factors have been accoimted for, 
the probable measurement error should be added .to accoxmt 
for differences in the target’s predicted bearing and range 
from the measurement system bearing and range values. 
Because of the uncertainty of these many factors, it is 
necessary to use wide bearing and range gates to insure 
that the target is not lost. The \mfortunate outcome of 
using wide bearing and range gates is the increased likeli- 
hood that two adjacent targets will enter each others 
bearing and range gates and become confused. 

The best solution to this dilemma is to allow the 
operator to select a target classification from his control 
station which would determine which of several predeter- 
mined bearing and range gate widths or width determining 
algorithms would be used with the designated target. As 
will be shown later, the designation of target maneuver- 
ability by the operator is not only useful in the measure- 
ment circuitry but also enhances the Kalman filter. 

The determining factor in the accuracy to which the 
bearing, range and time measurements can be represented 
in the microprocessor is the value represented by the 
least significant bit in the binary data word. If the 
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digitizing process is carried out with the digital zero 
equal to the analog zero it will introduce a mean error 
equal to one-half of the value represented by the least- 
significant bit and a standard deviation of error equal to 
3”^ times the value represented by the least significant 
bit. If the digital zero is offset by minus the previous 
mean error, then the mean error will be zero and the 

_ i 

standard deviation of error equal to 12~® times the value 
represented by the least significant bit. Obviously, the 
zero offset approach is preferred. 

The length of the binary data word used to represent 
the bearing, range, and time measurements is determined 
by the value represented by the least significant bit, 
and the maximum value which is to be represented. 

For the system under consideration, the least signifi- 
cant bit should represent the standard deviation of error 
of the target bearing with gyro error (O.I 258 d’egrees or 
0.002196 radians) and the target range ( 7 . 057 ^ yards), 
and the useful precision of the time measurement (+ 0.784 
msec), as a minimum. The maximum value of each measurement 
should be 2 ‘Tf radians, 60,000 yards and approximately 8.0 
secs. This would require 11.5 bits for the bearing measure- 
ment, 13.1 bits for the range measurement and 13*3 bits 
for the time measurement. If an 8-bit microprocessor is 
to be used, then data words of 8, 16, or 24 bit lengths 
are practical. In addition, the INTEL 8 O 8 O microprocessor 
has some instructions for manipulating l6-bit data words. 
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Thus 16-bit data words were considered optimum. This 
results in the following data format: 



Measurement 

Bearing 

Range 

Time 



Range 
0-2-tf rad 
0-60,000 yds 
0-8.0 sec 



MSB 
Tf rad 

30.000 yds 

4.0 sec 



LSB 

0.0959 mrad 

0.9155 yds 

0.1221 msec 



By changing the most significant bit to the next even power 
of two, the units can be rationalized as follows: 



Measurement 

Bearing 

Range 

Time 



Range 
0-2*Tr rad 

0-65,535 yds 

0 - 7«999 sec 



MSB 

4 . 0 rad 
32,768 yds 
4.0 sec 



LSB 

2"^^ rad 



1 . 0 yds 
2”^^ sec 



The resolution resulting from this format will result in 
less than 1.0 percent increase in the standard deviation 
of error as a result of digitizing the data. The remaining 
question is what resolution is possible with the available 
hardware . 

The bearing measurement could be conveniently made 
by combining one of the synchro-to-digital converter modules 
such as Transmagnetics series I 65 I, which has l4-bit 
resolution, and two l6-bit latchs , enabled by the video 
circuitry, to store the first and last bearings from 
which video was detected from a particular target on 
each scan. The microprocessor software could then compute 
the center bearing to the target. An alternate scheme 
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would "be to use a single latch to store the first hearing 
from which video was detected from a particular target on 
each scan and a counter to count the number of video pulses 
from the target on each scan. The microprocessor software 
could then use this data to estimate the center bearing. 

A third and far less accurate scheme would be to simply 
use the bearing of the first video pulse. The first 
scheme is recommended. 

The rainge measurement is essentially a time interval 
measurement. A range of 1 yard corresponds to 6.1 nsec. 
Thus a binary counter would have to be clocked at 163*9 
MHZ to count radar range in intervals of 1 yard. Ordinary 
Transistor-Transistor-Logic (TTL) counters are limited to 
clock frequencies of 30 MHZ and below, while high-speed TTL 
counters can operate as high as 50 MHZ Z~10_7. Schottky- 
Clamped TTL counters can operate reliably at clock fre- 
quencies of 100 MHZ Z~10_7. Emitter-Coupled-Logic (ECL) 
circuitry could perform at higher clock frequencies but 
the complexity and cost of the circuitry would increase 
considerably due to the unavailability of ECL Medium-Scale 
Integrated (MSI) semiconductor chips with the desired 
functions. A Shottky-C lamped TTL Counter, such as the 
74S197» is suggested, with a compromise reduction in the 
clock frequency to 81 . 96 IVIHZ resulting in a minimum range 
resolution of 2 yards. An alternate scheme would be to 
mix ECL circuitry with the Shottky counters. 

Another consideration in the range measurement counter 
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is the propagation of carry pulses through the counter. The 
7^3197 binary counter has a maximum propagation time of 13 
nsec through each stage. For the 15 stages required for a 
maximum range of 65 » 536 yards , a total propagation time of 
195 nsec would be required for a carry propagation to the 
highest significant bit. This is 16 times the clock period, 
which would lead to unstable data at the counter outputs. 

The easiest solution to this dilemma is to stop the clock 
input to the counter at the time the target video is 
detected and store the counter value after it has stabilized. 
This would necessitate a separate range coimter in the 
interface hardware for each target being tracked. 

The time measurement interface is comparatively straight 
forward. A l6-bit counter driven by a 8,192 HZ clock, with 
a storage latch for each target will suffice. As in the 
range counter, the carry propagation delay must be taken 
into account when storing the time measurement. Due to 
the slower clock frequency, the counter output will be 
stable 99 percent of the time. Simple combinational logic 
circuitry could be employed to synchronize the time measure- 
ment with the clock to avoid storing unstable counter 
outputs . 

Using currently available hardware, the radar interface 
can be accomplished with resolutions well within the 
radar measurement errors specified earlier in this section. 
The following data formats are recommended. 
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Measurement 


Ran^ 


No. of Bits 


Resolution 


Bearing 


0-2tt rad 


14 


•3835 mrad 


Range 


0 - 65,535 yards 


15 


2 . 0 yards 


Time 


0-8 sec 


16 


.1221 msec 


B . OWN SHIP 


COURSE AND SPEED 


INTERFACE 





The own ship course and speed interface serves two 
fimctions; (1) to smooth the variations in own ship's 
course and speed, and (2) to measure and store the smoothed 
course and speed to be associated with a particular target. 
There are two options of how to perform these functions, 
stemming from the fact that the microprocessor software 
uses own ship's course and speed data solely to compute 
own ship north/south and east/west velocities. 

The first approach would be to process the course and 
speed measurements independently. The course measurement 
could be smoothed either mechanically or with an over-damped 
servo feedback system. A synchro to digital angle converter 
could be employed to digitize the measurement. The speed 
measurement could be smoothed electrically and digitized 
using an analog to digital converter. 

The second approach would be to combine the course and 
speed measurements into north/south and east/west velo- 
cities. The course synchro signal could be converted to 
two analog signals corresponding to the sine and cosine of 
own ship's course using a synchro to DC sine/cosine converter, 
such as Transmagnetics model 655N V (series). These analog 
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signals could then he multiplied hy an analog speed signal 
with two DC analog multipliers. Low-pass filters could be 
used to smooth the two signals and then digital to analog 
converters used to digitize them. This approach would 
simplify the work of the microprocessor but additional 
errors would be generated in the analog multipliers . 

A compromise approach would be to use the first system 
but substitue a synchro to digital sine/cosine converter 
in place of the digital angle converter. With this tech- 
nique, the microprocessor would receive smoothed north/ 
south and east/west course component measurements and 
smoothed speed measurements . 

Assuming a standard deviation of error in the course 
measurement of 0.1 degree (l.?45 x 10~^ radians) prior 
to digitizing, a digital data format of 11.8 bits would be 
required for minimum resolution. Using a data format 
similar to that used for the target bearing (MSB = 2.0 
radians) and l4-bit resolution, the standard deviation 
of error will only increase 0.3 percent as a result of 
digitizing. Assuming a standard deviation of error in the 
speed measurement of 0.1 knot prior to digitizing, and a 
maximum speed of 100 knots, a digital data format of 10 
bits would be required. Using a 12-bit analog to digital 
converter or a 12-bit synchro to digital converter with the 
most significant bit representing 64 knots, the standard 
deviation of error will increase only 1.6 percent, as a 
result of digitizing. 
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In summary, the own ship course and speed interface 
could he accomplished using available hardware, with the 
following data fomats and measurement precision. 



Measurement 


Range 


No. of Bits 


Resolution 


Course 


0-2'rfrad 


14 


.3835 mrad 


Speed 


0-128 knots 


12 


.03125 knots 


C . OPERATOR/USER INTERFACE 







The operator/user interface serves two functions; (1) 
to provide operator control of the tracking system, sind 
(2) to display the output data to the users. 



The operator should have control over (1) which radar 
video is to be tracked, (2) the number of targets being 
tracked, and (3) the tracking parameters to be used with 
each target. 

It is proposed that the operator would control the 
system in the following manner. The measurement interface 
circuitry would be divided into tracking elements, each 
capable of handling one target. Each element would be 
composed of a bearing and range gate generation circuit, 
a range measurement circuit, and storage circuits for the 
bearing, range, time, own ship's course and own ship's 
speed data. 

The operator's control panel would have an off/on/start 
switch for each tracking element. The system would be 
turned on with all tracking elements turned off . To track 
a target, the operator would place his radar repeater's 
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bearing and range cursor on the target video to be tracked 
and momentarily place the tracking element switch in the 
start position when the radar sweep passed the target and 
then place the tracking element switch in the on position. 
This process would momentarily couple the radar repeater's 
bearing and range cursor video to the tracking elements 
bearing and range gate circuit and enable the tracker's 
measurement circuitry, thus taking a measurement of the 
target's position. If the measurement circuitry detected 
target video, an indicator light would be lit. To stop 
tracking a target, the operator would simply turn the 
tracking element off. A control data word (8 bits) would 
be periodically read by the microprocessor to advise it as 
to which tracking elements have updated measurement data 
available. The operator would have a second multi-position 
switch which would allow him to select the bearing and 
range gate widths and the variance of acceleration para- 
meters to be used with each target. 

The output data to be displayed would include (1) target 
designation, (2) target bearing, (3) target range, (4) tar- 
get course, and (5) target speed, as a minimum. For more . 
sophisticated systems the target's closest point of approach 
(CPA) could be computed and displayed as (1) CPA bearing, 

(2) CPA range and (3) time of CPA. The speed of the data 
output is critical to the capability of the system to track 
multiple targets. Direct memory access techniques are 
suggested as the only practical solution. 
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After the display data is extracted from the tracking 
system, several display techniques could be used. Tabular 
displays in the form of 7-segment LED modules or cathode 
ray tubes would facilitate displaying the data to several 
users at the operators station and also transmission of 
the data to remote stations, such as the bridge, CO's 
cabin, and weapons stations. One of the major advantages 
of an automatic tracking system is its capability to 
instantly and continuously disseminate tactical data to 
remote stations thus eliminating the need for passing 
this information on sound powered phones as is currently 
the practice on non-NTDS ships . 
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IV. KALMAN FILTER 



As stated earlier, the Kalman filter is a common optimal 
filtering technique for estimating the state of a linear 
system. "Among the presumed advantages of this type of 
data processor are that it minimizes the estimation error 
in a well defined statistical sense and that it utilizes 
all measurement data plus prior knowledge about the system. 
The corresponding potential disadvantages are its sensitiv- 
ity to erroneous a priori models and statistics, and the 
inherent computational burden." _/ 

A. DISCRETE KALMAN FILTER ALGORITHM 

The rotation of the surface search radar antenna 
results in periodic measurements of the target positions. 
Accordingly, the discrete form of the Kalman filter, rather 
than the continuous form, is appropriate. The discrete 
Kalman filter can be described by a set of matrix 
equations . 

Using the target state (Z^) representation presented in 
the description of the tracking problem, in Section I B, 
the matrix form of the measurement process is /~3_7s 



w = HZ + V 
n n n 



w 



n 






‘*5 



Xjjj(n) = Measured x Coordinate at Scan n 
y-j^(n) = Measured y Coordinate at Scam n 



H = 



10 0 0 
0 0 10 



= Target State Vector at Scan n 



V = 

n 



v^(n) 

Vy(n) 



V (n) = Random Measurement Error on x (n) 

X. in 

Vy(n) = Random Measurement Error on 
The estimated state vector at scan n (Z ) is 

= S ^n-1 

i - (See Tracking Problem Page 14) 

Z ^ = Optimal Estimate of Target State Vector 
at Scan n-1 



The optimal estimate of the target state at scan n (Z^) 



is given by Z~3_7: 



n 






, + K 

n n 


( w 
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■ H 


’XX 


V 




! /T 

XX ' 


®xy 
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‘yx 
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! /T 

yx ^ 


®yy 


/T 



n 



= Gain Matrix at Scan n 



The gain matrix (K^) is computed from /~3_7* 



rn rp -t 

K = P H^(HP + R ) ^ 
n n ' n n^ 
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P = Estimated Covariance of Estimation Errors 
^ Matrix Prior to Processing Scan n Measurements 




a/(n) 




o-/(n) 



= Covariance of Measurement 
Error Matrix at Scan n 



Following the processing of scan n measurements, an 
optimal estimate of the covariance of estimation error 
matrix for scan n can be computed Z”3_7* 

P = (I - K H)P 
n n n 

I = Identity Matrix 



From this, it is possible to compute the estimate of 
the covariance of estimation error matrix for scan n + 1 

rsy. 

^n+l " *n+l *n+l ^ n+1 ® „+i 

2 

Q ~ cr = Variance of Random Target Acceleration 

3 . 



The basis for the discrete Kalman filter in this form 
is the legitimacy of the target maneuver model. The 
variance of target acceleration must be equal and indepen- 
dent along the x and y axis, and from scan to scan. This 
is true for the targets of interest. Thus, this fora of 
the discrete Kalman filter is valid if the position 
measurements were made in the x and y coordinates. 

Using the polar coordinate system shown in Figure 4, 
the following transformations can be used with the radar's 
bearing (9) and range measurements (R). 
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FIGURE 4 COORDINATE SYSTEM 
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X = R sin 0 
m 

ym = R cos 0 

Using these transformations, the elements of the covariance 
of measurement error matrix are Z~ll_7: 



<7-2 

r 


sin^0 


+ R^ 


°0 


cos 0 




cos^0 


+ R^ 


^ 2 
^0 


. 2 
sin 



i sin 2 0 

2 

cr^ = Variance of Range Error 
2 

o'q = Variance of Bearing Error 

This modification completes the description of the discrete 
Kalman filter. 

In some applications, it is possible to precompute the 
values of the gain matrix and store them in the computers 
memory for subsequent use. This greatly simplifies the 
software algorithm. There are two reasons why this is not 
possible in this case. 

The first is that the elapsed time between measurements 
(T) is not constant. This is primarily due to the fact that 
the radar antenna rotates relative to own ship’s heading. 
Thus, changes in own ship's course can significantly 
increase or decrease the elapsed time between measurements 
(T). 

The second reason is that the covariances of measurement 
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error (cr^, (Ty, o^y) are a function of the hearing and range 
to the target. Consequently, the gain matrix must be com- 
puted on line for each target being tracked. A flow chart 
of these calculations is shown in Figure 5* 

The solution of the Kalman filter equations in matrix 
form is unnecessarily complex for the microprocessor. 
Therefore, an algebraic solution for the value of each 
element in the various matrices was formulated. Appendix 
A lists this algebraic form of the Kalman filter equations. 

B. KALMAN FILTER PARAMETERS 

The performance of the Kalman filter is determined by 

three parameters? (1) variance of bearing error (Oq ), 

2 

(2) variance of range error (cr'^ ), and (3) variance of 

2 

target acceleration {cr ). The first two parameters 
are functions of the radar, the ship's gyro, and the inter- 
face system.’ The third parameter is a function of the 
target's relative motion and the desired smoothness of the 
output data. 

The variance of bearing error is equal to the sum of 
the variance of the radar bearing error, (.0058 deg ), the 
variance of the bearing- interface digitizing error (.00016 
deg ) , and the variance of the ship's gyro error (0.01 deg ) . 
Thus, the total variance of bearing error is 0.016 deg or 
.0000049 rad^. 

The variance of range error is equal to the si^m of the 

2 2 

variance of the radar range error (50 yd or 42 m ) and 
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FIGURE 5 DISCRETE KALMAN FILTER FLOW DIAGRAM 




the variance of the range interface digitizing error 
2 2 

( 1.3 yd or 1.1 m). Thus the total variance of range 

2 2 

error is equal to 5i yd (43 m ) . 

The variance of target acceleration parameter is a 
compromise between the user's desired smoothness of the 
output data for non-maneuvering targets and the ability of 
the filter to provide the closest estimate of a target's 
position and velocity for a maneuvering target. 

A possible solution to this dilemma is to allow the 
operator to select the variance of acceleration parameter 
to be used with each target. This parameter, along with 
the bearing and range gate widths, could be selected by 
enabling one of several read-only-memories (ROMs) containing 
the desired parameters. The particular ROMs installed 
at any one time could be selected from a library of such 
ROMs by the Commanding Officer depending on the ship's 
mission. 

The determination of what values of variance of accelera- 
tion to be used with various targets are best determined by 
modeling Z~4_7. A computer model of the tracking problem 
was written and used to investigate the system performance 
with various values of variance of acceleration. The 
results of this investigation are presented in the figures 
discussed in Section IV E. 

C. KALMAN FILTER INITIAL CONDITIONS 

As illustrated by Figure 5» there are two basic loops in 
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the Kalman filter algorithm. They represent the recursive 
calculations of the optimal estimates of the covariance of 
estimation error matrix (P) and the target state matrix (Z). 
An initial condition must he established for each of these 
based on a priori information of the measurement process and 
the targets of interest. 

The initial condition of the target state matrix is 
equal to the expected value of the initial state. The 
expected values of the initial position elements, Xq and 
Yq, are the first position measurements; 

Xq = R(0) sin 9 ( 0 ) 

Yq = R(0) cos 9(0) 

The expected values of the initial velocity elements, X^ 
and Yq, are zero. Thus, the filtering of the target state 
does not begin until the second measurement. 

The initial condition of the covariance of estimation 
error matrix is equal to the covariance of error of the 
initial state estimate. 

Po = P ^^0 - ^0-7 - 7o-7 " 

E^(Xq-Xq)^^ 0 0 0 

0 E/"(Xq)2^ 0 0 

0 0 E^(Yq-Yq)2J7 0 

0 0 0 

The non-zero elements in this matrix are functions of 
the initial relative velocities which might be encoimtered, 
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and the measurement errors of the initial position. Assuming 

a uniform probability distribution of relative velocities 

from -50 yd/sec to +50 yd/sec (+ 46 m/sec), and an initial 

range of 50 f 000 yds ( 46,000 m), the variance of initial 

2 2 

position error is 12,000 yds (10,000 m ) and the variance 
of initial velocity error is 830 yd^/sec^ (700 m^/sec^) 
for both the X and Y coordinates . 

The primary influence of the initial conditions is 
the time required for the filter to achieve a steady-state 
solution. Large values lead to longer settling times for 
nonmaneuvering targets where small values lead to longer 
settling times for maneuvering targets. 

D. SYSTEM COMPUTER MODEL 

In order to evaluate the performance of the Kalman filter 
with various targets and parameters, a computer model of 
the system was derived. The computer model included the 
simulation of target and own ship motion, the radar measure- 
ment and interface process, the Kalman filter, and the 
computation of the filtered target course and speed. 
Additionally, the program included computation of measure- 
ment errors, target track errors, and associated statistics. 

The computer model simulated the target and own ship 
motion by digitally integrating the north/south and east/ 
west velocity components using an integration step of 0.01 
sec and double precision arithmetic. Target and own ship 
course changes were accomplished by a second-order nonlinear 
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control system model with the ordered course as the input, 
a limit on the maximum angular rate, and the actual course 
steered as the output. The control system model parameters 
and angular rate limit were adjusted to simulate typical 
course changes for the targets of interest. Speed changes 
were accomplished by a first-order linear control system 
model with the ordered speed as the input and the actual 
speed as the output. No attempt was made to duplicate 
actual ship maneuvers , but only to approximate the maneuvers 
typical of the targets of interest. 

The radar and interface systems were modeled by first 
comparing the rotation of the radar antenna to the bear- 
ing to each of the targets. Upon detection of a target, 
the bearing and range to the target were noted along with 
own ship's course and speed, and the time of detection. 
Random numbers representing the Gaussian distribution of 
measurement errors were added to the bearing and range 
measurement. The measurements were then truncated to 
simulate the digitizing process of the interface system. 

The Kalman filter was implemented through FORTRAN 
statements similar to the algebraic expressions listed in 
Appendix A. The only modifications were those related to 
improving the efficiency of executing the computations. 

The target's position, course, and speed were computed 
from the relative north/south and east/west positions 
and velocities provided by the Kalman filter, and the 
simulated measurement of ov/n ship's course and speed. 
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The computer model was programmed in FORTRAN IV, com- 
piled with a FORTRAN H Compiler, and executed on the Naval 
Postgraduate School's IBM 360/67 computer. 

E. KALMAN FILTER PERFORMANCE 

Using the computer model previously described, the per- 
formance of the Kalman filter was evaluated for two classes 
of targets, maneuvering and nonmaneuvering. In all cases, 
own ship was proceeding on a steady course ( 000 °) and at 
a fixed speed of 20 knots (37 km/h) from an initial position 
at the origin of the cartesian coordinate system. Eight 
targets were initially positioned and proceeding on the 
courses and at the speeds shown in Table 1 . A unique set 
of random numbers was used with each target for the various 
computer runs. This eliminated the influence of variations 
in the random number distributions on the statistics of 
the Kalman filter output. 

The following system parameters were used: 



RADAR PARAMETERS 



Pulse Repetition Rate 
Range Precision 



100.0 pps 

7.057^ yds 
(6.453 m) 

0.1258 deg 

17.0 rpm 



Bearing Precision 



Antenna Rotation Rate 



DIGITAL INTERFACE PRECISION PARAMETERS 



Bearing LSB 



Range LSB 



4.0 yds 
(3.6576 m) 

0.02197 deg 
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TARGET N/S POSITION E/W POSITION COURSE SPEED 

NO. ( vards/meters ) ( yards/me ters ) (degrees) ( knots /kms ) 
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TABLE I INITIAL POSITION, COURSE AND SPEED FOR SIMULATED TARGETS 







I 







Time LSB 



0.1221 msec 



Own Ship Course LSB 0.02197 deg 

Own Ship Speed ISB O .03125 knot 

( 0.05788 km/h) 

p 

The Kalman filter measurement error parameters (cri and 

2 

0 ^ ) were selected to match the computed measurement errors 
resulting from the radar and interface parameters. 

cf^ = 4.870 X 10"^ rad^ (0.01599 deg^) 

= 51.14 yd^ (42.76 m^) 

The Kalman filter initial conditions were the same as those 
discussed in Section C of this chapter. The variance of 
target acceleration (cr ) parameter was chosen as the inde- 
pendent variable. 

Both scenarios began with the acquisition process for 
each of the eight targets, which involved setting the 
initial conditions of the Kalman filter, and ran for 5 
minutes of simulated time. For the nonmaneuvering case, 
no changes were made throughout the run. For the maneuver- 
ing case, each of the targets executed a 45 deg course 
change beginning 120 sec after the start of the run. 

Targets 1, 4, 6 and 8 turned right, while targets 2, 3> 

5 and 7 turned left. These turns were executed with an 
angular rate of 3*0 deg/sec. Figure 6 shows an example 
of the plot developed during such a left turn. 

The actual measurement error statistics, based on 680 
measurement samples from the computer model, for the 
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6 SAMPLED GEOGRAPHIC TRACK FOR TARGET 3 
DURING 45 DEGREE COURSE CHANGE 
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nonmaneuvering and the maneuvering target scenarios are 
listed in Table II. There was good correlation between the 
measured variances of error and the computed variances . 
Additionally, the mean errors were approximately equal to 
minus one-half the values of the least significant bit 
used in the interface simulation. This is in agreement 
with the theoretical mean error generated by truncation. 

Each scenario was repeated for seventeen logarithmically 
distributed values of the Kalman filter variance of target 

o li o 

acceleration parameter from 0.0005 yd /sec (0.00042 m / 
i}, 2 / 4* / 2 / 

sec ) to 1 yd /sec (0.85 ^ /sec ). In all cases, the 
Kalman filter achieved a steady state condition in 60 sec, 
corresponding to 17 target position measurements. Steady 
state is defined as the condition where the Kalman filter 
gain matrix (K^) would be constant for targets with constant 
relative position and constant sample intervals. The 
bearing, range, course and speed of each target was com- 
puted from the Kalman filter's optimal estimate of target 
state and compared with the actual target bearing, range, 

t 

course and speed. Error averages were computed for the 
time interval from 60 sec to 300 sec. Figures 7, 8, 9 
and 10 show plots of the average root mean square of error 
for the target bearings, ranges, courses and speeds, 
respectively, as a function of the variance of target 
acceleration parameter. It should be noted that the 
error averages for the maneuvering target are a function of 
the elapsed time between maneuvers and the time spent in 
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TABLE II MEASUREMENT ERROR STATISTICS 
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FIGURE 7 AVERAGE ROOT MEAN SQUARE OF BEARING ERROR VS. 

VARIANCE OF TARGET ACCELERATION PARAMETER (o' 
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FIGURE 8 AVERAGE ROOT MEAN SQUARE OF RANGE ERROR VS. 

VARIANCE OF TARGET ACCELERATION PARAMETER (o' 
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FIGURE 9 AVERAGE ROOT MEAN SQUARE OF COURSE ERROR VS. 

VARIANCE OF TARGET ACCELERATION PARAMETER (a; 




t 

o 

o 



I 1 



-r 



m 



^ U 
<D O O 
<D U ^ 



1 ^ 

o 



o 

o 

o 



o 



65 



FIGURE 10 AVERAGE ROOT MEAN SQUARE OF SPEED ERROR VS. 

VARIANCE OF TARGET ACCELERATION PARAMETER (cr 



maneuver. Thus, the plots for the maneuvering target only 
apply for the particular scenario used. 

Several conclusions can be dravm from these figures. 

For nonmaneuvering targets, the variance of acceleration 

p oh, 

parameter (cr ) should be small (< .001 yd /sec (.0008 

3 . 

2 / 

m /sec )) if smooth course and speed values are desired. 

For maneuvering targets, there is an optiaml value of 
O' ^ (».10 yd^/sec^ (.08 m^/sec^)) below which the filter 

cL 

does not follow the target maneuvers well and above which 
the target data will be noisy regardless of the target 
maneuvers . There is a lower limit to how well the Kalman 
filter will track a maneuvering target even with an optimal 
value of O' . Finally, the bearing and range to a non- 

3 , 

maneuvering target will be improved by the Kalman filter, 
as shown by the standard deviation of measurement error 
line on Figures 7 and 8, over a wide range of o* values. 

3 . 
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V. MICROPROCESSOR PROGRAM 



A MCS-8 microcomputer system, using INTEL’S 8080B 
microprocessor, was used to implement the Kalman filter 
algorithm. The objective of implementing the algorithm 
was to determine limits to the time required to accomplish 
the Kalman filter computations and the accuracy of the 
results . 

A. PROGRAMMING LANGUAGE 

INTEL Corporation has developed the PL/M compiler, a 
derivative of PL/l, for use with their 8080 microprocessor. 
PL/M is a block structured language which speeds micro- 
processor program development by relieving the programmer 
of the tedious task of memory and register management 
It facilitates the use of variable names while still permit- 
ting the programmer the freedom to reference specific 
memory locations. The basic instruction set includes 
arithmetic and Boolean algebra operation on 8-bit and l6-bit 
variables and/or constants. Program control is accomplished 
with Procedure, Do, Go To, If, Call and Return statements. 
Logic functions such as comparison, shift and rotation are 
also included. Floating point and transcendental functions 
must be implemented by subroutines . PL/M permits the 
program control provided by assembly languages while 
taking a much more readable form. 
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B. FLOATING POINT ARITHMETIC 



Floating point arithmetic is required whenever the range 
of numbers represented by the program variables is large or 
unpredictable Such is the case with the Kalman 

filter algorithm. This is particularly true if the operator 
is given the freedom to alter the filter parameters . A 
disadvantage of using floating point arithmetic is the 
additional execution time required. 

A suitable set of floating point arithmetic procedures 
were not available. The procedures listed in Appendix B 
were derived from several existing procedures and the 
general guidelines detailed in Reference 12. 

A l6-bit mantissa, 7-bit exponent and 1-bit sign were 
chosen to represent the floating point numbers using three 
8-bit data words in the following format: 



7 0 7 0 7 0 



M 




M + 1 




1 M + 2 



Mantissa (l6-bits) ' Sign (1-bit) 

Exponent ( 7 -bits ) 

M - Memory Address Associated with the Floating Point 
Number 

The l6-bit mantissa was chosen because it is the largest 
binary mmiber that INTEL'S 8080 microprocessor can manipu- 
late without a considerable increase in the comjTLexity of 
the progi^Linming . Additionally, it is compatible with the 
input and output data. The mantissa is always left justi- 
fied with the binary point located to the left of the most 
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significant bit. Thus the mantissa is always equal to or 
greater than 0.5» and less than 1.0. 

The sign was located in the third word in order that the 
mantissa be a full l6-bits. A positive sign is represented 
by a binary 0 and a negative sign is represented by a 
binary 1 . 

The remaining 7-bits of the third word are used for the 
exponent. The value of the exponent is equal to 64 plus the 
power of two to be associated with the mantissa. 

The data and constants in the PL/M programs are written 
in this format with hexadecimal numbers. The following 
are examples of this form. 



DECIMAL 

0.5 

1.0 

- 1.0 

10.0 

64.0 

0.1 



HEXADECIMAL FLOATING POINT 
80H, OOH, 40H 
80H, OOH, 41H ^ 

80H, OOH, CIH 
AOH, OOH, 44H 
80H, OOH, 47H 
CCH, CDH, 3DH 



The largest magnitude number that can be represented 

is FFH, FFH, 7FH, which corresponds to 9-223 x 10^®. The 

smallest magnitude number that can be represented is 80H, 

— 20 

OOH, OOH which corresponds to 2.710 x 10 . Zero is 

represented by OOH, OOH, OOH. 

The typical execution time of the arithmetic and trig- 
onometric procedures listed in Appendix B, using INTEL'S 
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8080B microprocessor, are listed below. 



FLOATING POINT PROCEDURE MAXIMUM EXECUTION TIME (msec) 



Addition (ADD) 1,3 
Subtraction (SUB) 1.4 
Multiplication (MULT) 3.4 
Division (DIV) 4.1 
Comparison (COB/CPARE) 0.8 
Square Root (SQRT) 17*0 
Cos (TRIG 1) 18.0 
Sin (TRIG 2) 18.0 
Cos and Sin (TRIG 3) 30.0 
Arctangent (TRIG 4) 19.0 



The actual execution times of these procedures are a fionc- 
tion of the numbers involved. 

Overflow and underflow provisions are included in the 
arithmetic procedures. Overflow will result in a solution 
equal to tlie largest possible number with the appropriate 
sign. Underflow will result in a solution equal to zero. 

l/jC'jiHi u' 

The precision of the sine and cosine functions, using 
the procedure "TRIG”, is at worst + 0.000488 if the argu- 
ment is in the range 0 to Str. The precision of the arc- 
tangent function, using the procedure "TRIG", is equal 
to + 0.0610 mrad. 

The four arithmetic procedures are used by passing the 
variables' addresses to the procedure with a call statement. 
The following PL/M statements are examples. 
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FUNCTION 



pl/m statements 



X + Y = Z 
X - Y = Z 
A X B = C 
L f M = Q 



Call ADD ( .X, .Y, .Z) ; 
Call SUB ( .X, .Y, .Z) ; 
Call MULT ( .A, .B, .C) ; 
Call DIV ( .L, .M, .Q) ; 



Note: The dot symbol proceeding a variable name in PL/M 



indicates the address of the variable. Some 
variables may correspond to addresses of other 
variables, in which case the dot would be omitted. 



The "COMPARE" procedure returns a binary integer quantity 
which indicates the relationship of the arguments . 

pl/m statement 

N = COMPARE ( .A, .B) ; 

RESULT 

If A < B Then N = 0 

If A = B Then N = 1 

If A > B Then N = 2 



The argument for "SQRT" must be positive. 

FUNCTION PL/M STATEMENT 

Y Call SQRT ( .X, .Y) ; 

The procedure "TRIG" can be used to compute the sine 
and/or cosine, and arctangent by passing the addresses of 
the arguments plus a function indicator flag to the pro- 
cedure with a call statement. It requires less time to 
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compute sine and cosine of the same angle with one call 
rather than two. The following are examples of the multiple 
uses of the "TRIG" procedure. 



FUNCTION 
X = Cos (A) 

Y = Sin (A) 

X = Cos (A) and 

Y = Sin (A) 

A = Arc tan (Y/X) 



PL/M STATEMENTS 

Call TRIG ( .X, .W, .A, 1) ; 

Call TRIG (.W, .Y, .A, 2); 

Call TRIG (.X, .Y, .A, 3)j 
Call TRIG (.X, .Y, .A, 4); 



Note: W is a dummy variable used to fill the unused argu- 

ment location in the call statement when only cosine 
or sine is desired. 



C. TRACKING SYSTEM PROGRAM 

An experimental version of the tracking system software 
was written to determine the execution time of the Kalman 
filter algorithm and the error statistics of the resulting 
target track data. This program, written in PL/M, is 
listed in Appendix C. There are three functional parts 
to this program; (1) data input, (2) the Kalman filter, and 
(3) computation of output data. 

The data input portion of the program includes the 
simulated interface system output data, and the simulated 
control panel functions. The simulated interface system 
data was generated in the correct floating point format 
by the FORTRAN program used to model the system performance. 
The variables "BRG," "RNG, " and "TIM" correspond to the 
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measured target bearing, range and time of detection, 
respectively. The ovm ship north/south and east/west 
velocity are represented by the constants "NSVEL" and 
"EWVEL. " The Kalman filter parameters cr ^ and cr ^ 

SL w 3T 

are represented by "SIGA," "SIGTH," and "SIGR." The 
variables "MARKOLD" and "CONTROL" serve the functions of 
distinguishing which target data corresponds to existing 
tracks, and which control panel switches are turned on. 
These variables are associated with the variables "SIMFLAG 
( ) " and "DFLAG" which correspond to the interface unit 
from which the target data came and thus the target's track 
number. In an actual operating system the program would 
use INPUT statements or direct memory access to read in the 
target measurement and control panel switch data. The 
measurement data would have to be converted from fixed 
point to floating point format before processing. 

The Kalman filter statements are grouped into two long 
DO statements. The first is used to initialize the filter 
when a new target is first processed. The second group is 
used to process new data for existing target tracks. The 
sequence of statements and variable names parallel the 
algebraic expressions listed in Appendix A. Temporary 
variables were introduced where necessary to avoid repeti- 
tive operations and to store intermediate results. 

Using the Kalman filter optimal estimates of target 
position and relative velocity, the third portion of the 
program computes and displays the target bearing, range, 
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course and speed. 

Following the program debugging, a series of tests, 
using a full range of typical input data, were run to 
verify the correct computation of each variable. The 
FORTRAN model was used to compute the "correct" values for 
each variable using the Naval Postgraduate School IBM 360 
computer. It was found that the computational errors were 
confined to a range of significance sufficiently small to 
avoid instabilities in the filter. After processing several 
measurements, the Kalman filter estimates of target state 
computed with the microprocessor were within the same 
range of error as the IBM 360 solutions. The covariance 
of estimation error and Kalman gain matrix elements computed 
by the microprocessor converged rapidly on the values 
computed by the IBM 3^0 computer. 

As a check on the microprocessor's performance a series 
of eight test runs were made, each with a different target. 
The eight simulated targets were identical to the non- 
maneuvering targets used previously except the signs of 
the east/west initial positions for targets 3-8 were 
negative. This change was made because the execution time 
and precision of the "TRIG" procedure is worst for targets 
in the north-west quadrant. Nonmaneuvering targets were 
used because the available memory space limited the amount 
of data that could be loaded with the program. The test 
runs were restricted to I 30 sec tracks. A variance of 
acceleration of .001 yd /see (.00084 m /sec ) was selected 
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to insure that computational errors would not be over- 
shadowed by normal estimation error in the output. It 
was assxamed that these tests would represent worst case 
conditions . 

The test program occupied a total of 2E00^^ bytes of 
memory. This was divided as follows; 



Floating Point Procedures 


98^16 


Bytes 


Kalman Filter Algorithm 


12C7i5 


Bytes 


Bearing, Range, Course and 
Speed Computations 




Bytes 


Variable Storage 


383i6 


Bytes 



The error averages for the test data are listed in 
Table III. These error averages represent the performance 
of the microprocessor tracking program after the Kalman 
filter achieved a steady state condition, approximately 
50 sec after initialization. This resulted in a relatively 
small sample set of 20-21 samples for each target. The 
error statistics of the input data during this time are 
listed in Table IV. As compared to the measurement error 
statistics presented in Table II, there is less correlation 
with the theoretical measurement errors for the microprocessor 
test data. Therefore, the error averages listed in Table 
III should be considered typical only for a short period 
of operation rather than for long periods with all possible 
targets . 

An analysis of these results reveals the microprocessor's 
influence on the track data. By comparing the course and 
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TARGET BEARING RANGE COURSE SPEED 

NO. (degrees) (yards) (degrees) (knots) 

Mean RMS Mean RMS Mean RMS Mean RMS 
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TABLE III MICROPROCESSOR TRACKING PROGRAM ERROR AVERAGES 



TARGET BEARING RANGE 

NO. (degrees) (yards) 

Mean Std. Dev. Mean Std. Dev 
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TRACKING PROGRAM INPUT DATA 



:rn;;ir 




speed data for targets 3i 5 and 7 with the data for targets 
4, 6 and 8 it can be seen that for targets with similar 
course and speed the errors increase with range. Addition- 
ally, errors are larger for targets with small relative 
velocities. The combination of long range and small rela- 
tive velocities can be expected to result in larger errors 
due to the limitations of the floating point arithmetic to 
handle large numbers with small differences. 

The final parameter of interest is the execution time 
of the microprocessor program. The INTEL 8080B micro- 
processor was able to execute 85 consecutive iterations of 
the Kalman filter algorithm in 37*5 secs, resulting in an 
average execution time of 0.44 secs. The computation of 
target course and speed required a worst case execution 
time of 54 msec. It would not be necessary to compute 
the target course and speed for every target on each antenna 
scan. The additional processing required to output the 
course and speed data would best be accomplished by a 
separate microprocessor dedicated to each display device. 

In general, a network of microprocessors could be used to 
share the task and thus increase the number of targets 
which could be tracked. 
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VI. CONCLUSION 



The tedious, error-prone task of tracking surface search 
radar targets manually, as is common onboard naval ships 
not equipped with NTDS, can be eliminated with the use of 
large scale integrated microprocessors. 

The AN/SPS-10 Surface Search Radar, installed on most 
U.S. Navy ships, could provide target position measurements 
with sufficient accuracy to perform automatic tracking of 
surface radar targets. Using leading edge detection and 
currently available hardware it would be possible to 
digitally measure the bearing and range to a target with 
a standard deviation of error of 0.126 deg and 7*15 yards 
( 6 . 5 ^ m) respectively. 

Using a 4-state Kalman filter and appropriate conver- 
sions for polar coordinate measurements, it is possible to 
obtain optimal estimates of target north/south and east/ 
west position and velocity. From this data and knowledge 
of own ship's motion, target course and speed can be com- 
puted far more accurately than by manual tracking tech- 
niques . 

A set of floating point arithmetic procedures was 
developed for INTEL'S 8080 microprocessor based on a 
floating point number format with a l6-bit mantissa and 
a 7-bit exponent. These procedures maintain sufficient 
accuracy to reliably compute the Kalman filter variables 



79 



while limiting the execution time sufficiently to pennit 
complex functions, like the Kalman filter, to he implemented. 

The Kalman filter, and course and speed calculations were 
successfully implemented using the microprocessor develop- 
ment system at the Naval Postgraduate School. ^ The average 
execution time of the Kalman filter was 0.44 secs while 
the course and speed calculations required 54 msec. This 
was accomplished with computed course and speed errors for 
nonmaneuvering targets typically less than 2.0 deg and 0.5 
knots. With the 3*5 sec sweep of the AN/SPS-10, a single 
system could continuously provide course and speed data 
for seven targets. This could he increased hy using faster 
microprocessors currently entering the market. 

It is recommended that future work in this area he 
devoted to (1) improvements in the floating-point software, 
( 2 ) development of the interface hardware, and (3) improve- 
ments to the Kalman filter algorithm. The floating point 
software could he improved hy reducing the execution time 
of the arithmetic procedures and increasing the precision 
of the transcendental functions. The Kalman filter algorithm 
could he improved hy minimizing the effects of own ship's 
motion, and reducing the execution time. 

A microprocessor-hased tracking system has the potential 
to provide any ship with a digital automatic surface track- 
ing capability, without the expense of NTDS size equipment. 

If implemented, it could reduce the manning of underway 
watches while improving the quality of the tactical 
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information available to the Officer of the Deck. There 
are additional applications of similar systems for such 
functions as sonar tracking and navigation. 
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P(n) = i P(n-1)«^ +PQr 
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APPENDIX A ALGEBRAIC FORM OF KALMAN FILTER 
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K(n) = P(n) H^(H P(n) + R(n) )“ 
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P(n) = (I - K(n) H) P(n) 
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APPENDIX B FLOATING POINT SOFTWARE 
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FLCATING POINT DIVICE ROUINE 
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